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The fluxes of atmospheric muons and neutrinos are calculated by a three dimensional Monte Carlo 
simulation with the air shower code CORSIKA using the hadronic interaction models DPMJET, 
VENUS, GHEISHA, and UrQMD. For the simulation of low energy primary particles the original 
CORSIKA has been extended by a parametrization of the solar modulation and a microscopic cal- 
culation of the directional dependence of the geomagnetic cut-off functions. An accurate description 
for the geography of the Earth has been included by a digital elevation model, tables for the local 
magnetic field in the atmosphere, and various atmospheric models for different geographic latitudes 
and annual seasons. CORSIKA is used to calculate atmospheric muon fluxes for different locations 
and the neutrino fluxes for Kamioka. The results of CORSIKA for the muon fluxes are verified 
by an extensive comparison with recent measurements. The obtained neutrino fluxes are compared 
with other calculations and the influence of the hadronic interaction model, the geomagnetic cut-off 
and the local magnetic field on the neutrino fluxes is investigated. 

PACS numbers: 95.85.Ry 96.40.Tv 



I. INTRODUCTION 

Atmospheric neutrinos are produced by the interaction 
of the primary cosmic radiation with the Earth's atmo- 
sphere. They result mainly from the decay of the charged 
pions and muons: 



7T — ► fl + 

<-> e + + v e + vjl 



(1) 



and to about 10 % from similar reaction chains for kaons. 
A simple balancing of the different neutrino species in- 
volved results in the following approximate relations be- 
tween the numbers of neutrinos: 
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1, and 



(2) 



A more detailed calculation leads to an energy depen- 
dence of all the ratios. Especially the ratio of muon 
neutrinos to electron neutrinos strongly depends on the 
energy, because the number of muons reaching sea-level 
before decaying increases with the energy. 

A precise simulation of atmospheric neutrino fluxes is 
of essential interest for the interpretation of the so-called 
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atmospheric neutrino anomaly, i.e. the observation of 
several neutrino detectors [Ell, El II, 01 Ll that the ratio of 
muon neutrinos to electron neutrinos in the atmosphere 
differs approximately by a factor of 2 from the theoretical 
predictions. The flux of electron neutrinos seems to agree 
relatively well with the expectation, and the anomaly 
results mainly from a lack of muon neutrinos. 

In addition, the anomaly displays a pronounced de- 
pendence on the angle of incidence. The highest deficit 
is measured for neutrinos traveling through the Earth 
for entering the detectors in upward direction, while for 
downward going neutrinos agreement with the theory is 
found. This directional dependence of the anomaly is 
commonly interpreted in terms of neutrino oscillations. 

Due to the enormous size of the detector, the re- 
sults obtained by the Super-Kamiokande experiment 
near Kamioka, Japan are statistically most significant 
and allow a most detailed exploration of the anomaly. As 
only detector so far, Super-Kamiokande could establish 
a pronounced East- West-effect in the neutrino flux orig- 
inating from the influence of the Earth's magnetic field 
on the trajectories of the charged primary and secondary 
cosmic ray particles 7]. 

The fluxes of atmospheric neutrinos have been calcu- 
lated with various theoretical approaches invoking differ- 
ent hadronic interaction models. Detailed calculations 
have been done by Barr, Gaisser, and Stanev (BGS) 0, 
E0; Bugaev and Naumov (BN) pj|; Honda, Kasa- 
hara, Hidaka, and Midorikawa (HKHM) 0, [T3|; Lee, 
Bludman, and Koh (LBK) HI 113, Tserkovnyak, Ko- 
mar, Nally, and Waltham (TKNW) [SJ El; Bat- 
tistoni, Ferrari, Lipari, Montaruli, Sala, and Rancati 
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TABLE I: Features of the different models applied in the calculation of atmospheric neutrino fluxes. The following abbreviations 
are used in the table: 1GRF stands for the International Geomagnetic Reference Field 0, WMM for the World Magnetic Field 
Model |glllO|. and USSA for the US Standard Atmosphere [Tl|. The terms used in the table are explained in Sec. ITU 



BGS BN HKHM LBK TKNW BFLMSR HKKM Ply 



Hadronic interaction TARGET semi- FRITIOF/ TARGET GEANT FLUKA FRITIOF/NUCRIN GEANT 
model analytical NUCRIN DPMJET III 

(parametrized) 



Dimensions 


1 


1 


1 


3 


3 


3 


3 


3 


Directional dependence 
of geomagnetic cut-off 


dipole 


dipole- 
like 


IGRF 


dipole 


IGRF 


IGRF 


dipole 


WMM 


Penumbra of cut-off 


no 


no 


no 


no 


yes 


yes 


no 


yes 


Local magnetic field 


no 


no 


no 


no 


yes 


no 


yes 


yes 


Energy loss by ionization 


yes 


yes 


yes 


yes 


yes 


yes 


yes 


yes 


Multiple scattering of muons 


no 


no 


no 


no 


yes 


yes 


yes 


yes 


Atmospheric model 


USSA 


Dorman 
model 


USSA 


? 


USSA 


USSA 


USSA 


USSA 


Elevation model of the 


no 


no 


no 


no 


no 


no 


no 


no 



Earth 



(BFLMSR) I2II2II2I: Honda. Kajita, Kasahara, and 
Midorikawa (HKKM)j2E E3 ; and Plyaskin (Ply) [23. A 
recent review of the calculations of atmospheric neutrinos 
can be found in Ref. psj . 

The calculation of BGS is a one dimensional Monte 
Carlo simulation made in two steps. First, cascades 
for different primary energies and zenith angles are sim- 
ulated, and subsequently, the energy depending yields 
of the secondary particles are weighted by the primary 
spectrum and the geomagnetic cut-off characteristics for 
the detector location. The hadronic interactions are de- 
scribed with TARGET [2^] , a parametrization of accel- 
erator data with special emphasis on energies around 
20GeV. 

The BN calculation is based on a one dimensional semi- 
classical integration of the atmospheric cascade equations 
in straight forward approximation over the primary spec- 
trum. The hadronic interaction is described by an an- 
alytical parametrization of double differential inclusive 
cross-sections based on a compilation of accelerator data. 
This approach neglects many details on the nature of the 
hadronic interaction [a]. 

The HKHM calculation is made by using the air 
shower simulation code COSMOS [3(J in a one dimen- 



[a] Note added in proof: After acceptance for publication we got 
knowledge about new more detailed numerical calculations of 
Naumov et al. I hep-ph/0201310 and ref erences therein) which 
essentially confirm the results of Ref. Il . 



sional Monte Carlo simulation. For energies above 5 GeV 
the hadronic interaction is described in the frame of 
FRITIOF version 1.6 [3l| with JETSET 6.3 [H At 
lower energies NUCRIN is used |33J. 

The model applied in the LBK calculations is the 
model of BGS, but extending the calculation on three 
dimensions. The same primary spectrum has been used, 
too. The calculation intends to study the influence of 
the transversal momenta in the different reactions on the 
neutrino fluxes. 

The three dimensional calculation of TKNW is based 
on the GEANT 3.21 detector simulation tool and its 
various models for the hadronic interaction, CALOR m, 
HIH, FLUKA92 HIIH, and GHEISHA HJ. 

Both, the LBK and the TKNW calculation failed to 
discover a major enhancement of the neutrino fluxes near 
the horizon, which was predicted for the first time in the 
three dimensional simulation of BFLMSR. In the mean- 
time, the TKNW group revised its model, and finds also 
an enhancement at the horizon |2lj . 

In the calculations of BFLMSR the FLUKA98 and 
FLUKA2000 are used as models for the hadronic 

interaction. These versions of FLUKA are quite different 
from the FLUKA92 version integrated in the GEANT 
package and being used in the TKNW calculation. 

HKKM extended the calculation of HKHM to three di- 
mensions. Additionally, the interaction models of COS- 
MOS can be replaced now by a parametrized version of 
DPMJET III g3|. Meaning that instead of interfacing 
DPMJET to COSMOS, DPMJET is run at fixed energies 
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and the yields of secondary particles are parametrized. 
This enhances very much the calculation speed, but sub- 
tile details of the interaction model might be lost in this 
approach. The published results being used in this paper 
for comparisons are based on the hadronic interaction 
models of the original COSMOS. 

The calculation of Plyaskin is based also on the 
GEANT detector simulation package, but only the 
GHEISHA model is used for the simulation of the 
hadronic interaction. The atmosphere is sampled in lay- 
ers with constant density of 1km thickness. 

The major differences of the various neutrino calcu- 
lations in handling certain physical effects and the geo- 
graphical details of the Earth are compiled in Tab. [I] 

In this communication a full three dimensional sim- 
ulation procedure for atmospheric muon and neutrino 
fluxes is presented using the standard air shower sim- 
ulation code CORSIKA [iif. In contrast to the previ- 
ous calculations which assumed the Earth for instance 
as mathematical sphere, the actual attempt includes a 
complete description of the geographical parameters of 
the Earth. For this purpose the CORSIKA 6.0 code is 
extended by a precise calculation of the geomagnetic cut- 
off, a parametrization of the solar modulation, a digital 
elevation model of the Earth, tables for the local mag- 
netic field in the atmosphere, and various atmospheric 
models for different climatic zones and annual seasons. 

As evident from Equs.^andElthe correlation between 
neutrinos and muons is very direct, especially the charge 
ratio of muons reflects the ratio of electron neutrinos to 
electron antineutrinos. Thus, the calculation results of 
CORSIKA are verified by the simulation of atmospheric 
muons and their comparison with recent measurements. 

The procedure of simulation is demonstrated by a 
detailed calculation of atmospheric neutrino fluxes in 
Kamioka. Using the versatility of the CORSIKA pro- 
gram to cooperate with different models for the simula- 
tion of the hadronic interaction, special emphasis is put 
on the question, how various formulations of the hadronic 
interaction influence the fluxes of atmospheric neutrinos. 
It will be shown that uncertainties in the description of 
the hadronic interaction are the main error source for the 
calculation of atmospheric neutrino fluxes. 

Furthermore, the effects of the geomagnetic cut-off 
modulating the primary flux and of the local magnetic 
field which deflects the charged shower particles on their 
way through the atmosphere are studied in detail. Re- 
peating the calculations, setting once the geomagnetic 
cut-off and once the local magnetic field zero, allows to 
disentangle the individual influences. 

It will be proven, that the local magnetic field is by far 
not negligible and leads to an increase of the ratio of elec- 
tron neutrinos to electron antineutrinos. Furthermore, 
it modulates the azimuthal dependence of the neutrino 
fluxes and causes an East- West-effect, clearly visible for 
sites with a low geomagnetic cut-off. 



II. THE SIMULATION TOOL CORSIKA AND 
ITS EXTENSIONS FOR THE SIMULATION OF 
LOW ENERGY ATMOSPHERIC PARTICLES 

A. The air shower simulation program CORSIKA 

The simulation tool CORSIKA has been originally de- 
signed for the four dimensional simulation of extensive 
air showers with primary energies around 10 15 eV. The 
particle transport includes the particle ranges defined by 
the life time of the particle and its cross-section with air. 
The density profile of the atmosphere is handled as con- 
tinuous function, thus not sampled in layers of constant 
density. 

Ionization losses, multiple scattering, and the deflec- 
tion in the local magnetic field are considered. The de- 
cay of particles is simulated in exact kinematics, and the 
muon polarization is taken into account. 

In contrast to other air shower simulations tools, COR- 
SIKA offers alternatively six different models for the de- 
scription of the high energy hadronic interaction and 
three different models for the description of the low en- 
ergy hadronic interaction. The threshold between the 
high and low energy models is set by default to Ehab — 
80GeV/n. 

Due to the steep spectrum of primary cosmic rays, 
only some 10% of the neutrinos detected in the Super- 
Kamiokande experiment originate from primary particles 
with energies higher than 80 GeV/n and the quality of the 
simulated neutrino fluxes mainly depends on the trusti- 
ness in the models describing the low energy hadronic 
interaction. 

Nevertheless, the extent to which different high energy 
interaction models are able to reproduce experimental 
muon data has been investigated in a previous paper |45| . 
It has been shown that DPMJET II.5 [M EH 111 and 
VENUS 4.125 agree best with muon data, while 
QGSJET [13 and SIBYLL 1.6 [U E3 are not repro- 
ducing well the charge ratio of muons above 80GeV/n. 
Thus, in this paper only DPMJET and VENUS are used 
for the simulation of the high energy hadronic interac- 
tion. 

For the simulation of the low energy hadronic inter- 
action GHEISHA H3 and UrQMD 1.1 [H 13 are ap- 
plied. Additionally DPMJET includes some extensions 
which allow also the simulation of the hadronic interac- 
tion down to energies of 1 GeV. In this case UrQMD is 
used for the simulation of hadronic interactions with en- 
ergies below. The total number of muons and neutrinos 
resulting from hadronic interactions below 1 GeV is very 
small, thus UrQMD plays more the role of a technical 
fallback, to prevent the program from crashing. No real 
influence of UrQMD is noticeable in the physical results 
in this case. 

Low energy reactions are handled very similarly in 
DPMJET II.5 and DPMJET III, so that the results ob- 
tained with both versions in the energy range relevant for 
the atmospheric neutrino anomaly should be fully com- 
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parable. 

Fluxes calculated by CORSIKA have statistical errors, 
caused by the limited number of particles calculated in 
the Monte Carlo simulation, and various systematic er- 
rors. It can be assumed that the main sources of sys- 
tematic errors result from the primary spectrum and the 
hadronic interaction models. Errors due to the particle 
tracking or particle decay are hardly to be quantified but 
they should be negligible compared to the other error 
sources. All errors given in the further results are purely 
statistical. 



B. The fluxes of primary cosmic particles 

A major uncertainty in the early calculations of at- 
mospheric particle fluxes stems from the absolute pri- 
mary particle fluxes. These are measured by satellite 
or balloon borne experiments, operating above or at the 
limit of the Earth's atmosphere. In Fig. the results 
of recent experiments are compiled. The balloon experi- 
ment MASS has been operated in Fort Sumner, New 
Mexico, where the vertical geomagnetic cut-off rigidity 
is 4.2 GV, explaining the missing flux below the cut-off. 
The balloon experiments BESS CAPRICE (57|, and 
IMAX [58| were launched in Lynn Lake, Canada, near 
the geomagnetic pole with a very low cut-off rigidity of 
about 0.5 GV. The Space Shuttle mission of the AMS 
prototype |5^, |6(| collected data over a large range of 
cut-off rigidities ranging from the maximum rigidity at 
the geomagnetic equator down to vertical cut-off rigidi- 
ties less than 0.2 GV, corresponding to proton momenta 
well below the pion production threshold. 

At low energies, especially below 10 GeV, solar modu- 
lation becomes important and introduces a further, time 
dependent source of differences between the data. But 
the experiments differ also at higher energies. The re- 
sults of AMS and BESS agree perfectly within experi- 
mental errors while all other experiments report fluxes, 
being mostly about 15-20% lower. The differences be- 
tween the experiments are not constant in energy and 
can therefore not be explained by a simple offset in the 
energy calibration. For instance the MASS and IMAX 
results agree at higher energies with the AMS data while 
the CAPRICE results match at lower energies. 

AMS and BESS detectors have been calibrated at ac- 
celerator beams of protons (BESS, AMS), He and C nu- 
clei (AMS). This ensures that the performance of the 
detectors and the analyzing procedure were thoroughly 
understood, giving high evidence that the higher primary 
proton fluxes reported by AMS and BESS are the better 
ones. 

The results for primary helium nuclei show similar dif- 
ferences between the experiments like the results for pri- 
mary protons. Again AMS and BESS report higher fluxes 
than the other experiments, but the results of AMS and 
BESS do not agree completely. The cross calibration 
with light ions in case of AMS is a strong argument for 
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FIG. 1: The fluxes of primary protons (a) and helium nuclei 
(b) as measured by recent balloon and satellite borne exper- 
iments. In order to enhance the differences in the region of 
interest, the fluxes are multiplied by E 2 ' 5 . 

the correctness of the AMS data. 

Different from the primary flux parametrization, pro- 
posed recently by Gaisser et al. |6l| . where the helium 
flux is obtained by a combined fit of the AMS and BESS 
results, the calculations in this paper are based on the 
AMS results, only. The primary particle generator in 
CORSIKA uses power laws extracted from the higher 
energy data of AMS including the solar modulation and 
the geomagnetic cut-off as described in the next sections. 

The bulk of primary particles producing neutrinos with 
energies being detected in Super-Kamiokande is covered 
by the momentum acceptance of AMS. In order to avoid 
any artificial cut for higher energies, the power laws have 
been just extrapolated up to the knee region. As our 
knowledge of the cosmic radiation at higher energies is 
rather poor, this assumption is still in fair agreement with 
the measurements. 



C. The description of the solar modulation 

The sun emits a magnetized plasma with a velocity of 
100 - 200km/s 62] . To reach the Earth, galactic cosmic 
rays have to diffuse into the inner heliosphcre against the 
outward flow of the turbulent solar wind, a process know 
as solar modulation. Depending on the solar activity the 
lowest energy cosmic particles reach the Earth with a 
variable flux. 

For most places on Earth the geomagnetic cut-off alters 
the primary particle fluxes stronger than the influence of 
the solar modulation. Therefore the geomagnetic cut- 
off must be simulated in a detailed microscopic calcula- 
tion as described in Sec. Ill Dl while the solar modulation 
can be handled by the parametrization of Gleeson and 
Axford |63(. This parametrization is based on a spheri- 
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FIG. 2: The fluxes of primary protons and helium nuclei 
as obtained by the primary particle generator of CORSIKA 
including the solar modulation but no geomagnetic cut-off 
compared with the results of AMS. In order to pronounce the 
differences in the region of interest, the fluxes are multiplied 
by E 2 . 



cal symmetrical model in which the differential intensity 
J(r, E, t) for the total energy E in the distance r from 
the sun for the time t is given by 
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J( ^' t) = FlR Jlo°,E + m , (3) 



with Eq being the rest mass and tp(t) is a free, time de- 
pending parameter which can be interpreted as the en- 
ergy loss of a primary particle during its approach to the 
Earth. 

In principle, tp(t) can be deduced within theoretical 
models from the solar activity. Nevertheless for the cal- 
culation of the neutrino fluxes in Kamioka, ^(t) can be 
assumed as constant in time. The flight of AMS took 
place roughly in the middle of the data taking period 
of Super-Kamiokande, thus the values of ip for primary 
protons and helium nuclei are obtained directly by a fit 
of the function in Equ. [3] to the low energy part of the 
spectra measured by the AMS experiment. The resulting 
absolute primary particle spectra without considering the 
geomagnetic cut-off used for the primary particle genera- 
tor of CORSIKA are shown in Fig. EI The overall agree- 
ment is quite good, but a systematical deviation around 
10 GeV indicates that the used parametrization is not the 
best possible. Nevertheless the deviation remains mostly 
within the experimental errors and the highest discrep- 
ancy for a single point is found as 6%. The additional 
error caused by this in the atmospheric particle fluxes is 
quite small. 



FIG. 3: The mean geomagnetic cut-off and the width of 
the penumbra region for Kamioka, Japan. The width of the 
penumbra region is defined by the rigidity difference of the 
lowest momentum of an antiproton escaping to outer space 
and the highest momentum of an antiproton being trapped. 
The counting of the azimuth angle (f>, here and in all further 
plots, follows the convention used by the Super-Kamiokande 
detector, = 0° means looking to the South (the particle 
travels to the North), (f> = 90° to the East (the particle travels 
to the West), etc. The North direction here is defined as the 
geographical one. The angle between the geomagnetic and 
geographic North direction in Kamioka is -7.59°. 



D. The simulation of the geomagnetic cut-off 

The Earth's magnetic field has nearly the shape of a 
dipole field. The field is strong enough to deflect charged 
primary particles on their way to the Earth's surface. 
While near the geomagnetic poles particles with very low 
momenta can penetrate to the Earth's surface, protons 
with up to 60 GeV, impinging horizontally near the geo- 
magnetic equator are reflected back to space. 

The calculation of the geomagnetic cut-off is done in 
a Monte Carlo simulation of the possible particle tra- 
jectories in the so-called back-tracking method. Instead 
of tracking primary protons from outer space to the 
Earth's surface, antiprotons from the surface are retraced 
to outer space. This method has the advantage that it 
allows to calculate in a straight forward way a table of 
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allowed and forbidden trajectories. The entries in the ta- 
ble depend on the location on Earth, the arrival direction 
and the particle momentum. 

In detail the particle tracking starts at 112.83 km, the 
top of atmosphere as defined in CORSIKA. The influence 
of the local magnetic field in the atmosphere, including 
the deflection of charged shower particles is handled lat- 
eron by CORSIKA using the approximation of a homoge- 
nous field. 

The particle tracking is based on GEANT 3.21 ^ 
and the magnetic field is described by the International 
Geomagnetic Reference Field 8] for the year 2000. For 
the downward going particle fluxes the location where 
the primary particle enters the atmosphere is bound to 
the nearer surrounding of the experiment and the arrival 
direction is sampled in cells of a solid angle of 250 fisr. 
For upward going neutrinos the geomagnetic cut-off is 
calculated for 1655 locations, distributed nearly equidis- 
tant over the Earth's surface, and the angle of incidence 
for each location is sampled in cells of 48 msr. 

Instead of calculating a sharp cut-off, functions in mo- 
mentum steps of 0.2GeV/c up to a maximum momen- 
tum of 64GeV/c are evaluated. This procedure accounts 
for the penumbra region of the cut-off, i.e. the chaotic 
change from open and closed trajectories which can be 
observed in irregular magnetic fields, as in case of the 
geomagnetic field. 

As an example for the results obtained for a fixed de- 
tector location, the mean geomagnetic cut-off for par- 
ticles entering the atmosphere in Kamioka is shown in 
Fig. |31 Local irregularities of the magnetic field over 
Japan cause a remarkable strong deviation from the reg- 
ular shape expected for a magnetic dipole field. Assum- 
ing highly accurate Monte Carlo simulations and highly 
accurate measurements, this feature should be reflected 
in the zenithal and azimuthal dependence of the particle 
intensities in Kamioka. 

Kamioka has a very extended penumbra region which 
exceeds a width of 4GV in some particular directions. 
Details about the simulation of the geomagnetic cut-off 
and plots for other locations on the Earth may be found 
in Ref. ||. 

A check of the primary particle generator in CORSIKA 
with its assumptions for the solar modulation and the 
geomagnetic cut-off can be done by the recent results of 
the AMS-prototype mission j6^. Due to the inclination 
of 51.7° of the shuttle orbit, the space craft passes geo- 
magnetic latitudes from to more than 1 rad. 

The experimental spectra of downward going protons 
and helium nuclei can be compared rather directly with 
the results of the primary particle generator. Only a cor- 
rection for the altitude dependence of the geomagnetic 
cut-off has to be applied. The cut-off has generally the 
highest value at the surface of the Earth, decreases with 
the altitude and vanishes when leaving the Earth's mag- 
netosphere. The mean difference in the cut-off between 
the top of atmosphere as assumed in CORSIKA and the 
orbit of the space shuttle is evaluated by a dedicated 
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FIG. 4: Comparison of the AMS-results on downward going 
primary protons for different intervals of the geomagnetic lat- 
itude with spectra produced by the primary particle genera- 
tor of CORSIKA including the simulation of the geomagnetic 
cut-off and the parametrization for the solar modulation. 



GEANT simulation and has a value of about 10%. 

The spectra of primary protons for different regions of 
the geomagnetic latitude together with the spectra pro- 
duced by the primary particle generator of CORSIKA 
are shown in Fig.0] The agreement between experiment 
and simulation is very good and the systematic decrease 
of the geomagnetic cut-off with the geomagnetic latitude 
is reproduced nicely. Only the spectrum for geomagnetic 
latitudes 0.9 < 9 mag < 1 shows a noticeable difference, 
which has to be attributed to the low absolute value of 
the cut-off which becomes comparable to the momentum 
steps used in the simulation of the cut-off functions. This 
disagreement has no significance for the calculation of at- 
mospheric muon or neutrino fluxes, because the primary 
energies are already near or below the pion production 
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threshold. The results obtained for primary helium nu- 
clei have a similar quality. 

Particles stored for longer times in the geomagnetic 
field, the so-called albedo or sub-threshold particles are 
not considered in the pres ent calculations. It has been 
demonstrated in Ref. J66J that they contribute to the 
atmospheric particle flux only negligibly. 



E. The geography of the Earth in CORSIKA 

The geography of the Earth plays a certain role in the 
simulation of atmospheric particle fluxes, because the ap- 
parent thickness of the atmosphere is altered by the dif- 
ferent elevation of the terrain over sea-level and various 
climatic conditions on the density structure of the atmo- 
sphere. Also the local geomagnetic field, bending charged 
secondary particles in the atmosphere has quite a differ- 
ent strength for locations near the geomagnetic poles and 
the equatorial regions. For the geomagnetic poles the ab- 
solute field is found to be 64.6 /iT, while the strength at 
the geomagnetic equator is only 21.7 /iT. 

Due to the vicinity of the place where the primary 
particle enters the atmosphere and the place of detection 
in the simulation of vertical downward going neutrinos, 
the geographic data are assumed as constant in the cor- 
responding calculations. For the simulation of inclined 
particles, the distance in the locations may reach already 
1200 km, and in case of upward going neutrinos the ori- 
gin of the primary cosmic particles is distributed over the 
entire Earth. 

Therefore, the local geomagnetic field is tabulated 
on basis of the International Geomagnetic Reference 
Field 8] by a table containing the field parameters for 
64800 locations distributed over the Earth's surface. The 
elevation over sea-level is described in a table of equiv- 
alent resolution of data published by the US National 
Geophysical Data Center in Ref. |67| . 

The atmospheric profiles observed in tropical and po- 
lar regions show considerable differences. The non-tropic 
atmospheres are subject to additional variations with the 
annual seasons. The extended CORSIKA code accounts 
for these effects by 7 atmospheric models 68j. The cor- 
responding density distributions are plotted in Fig.|3J As 
expected, the largest differences appear between the po- 
lar winter and summer. The seasonal variations become 
less important and vanish as the climatic zone approaches 
the equator. 



F. The settings and the way of simulation in 
CORSIKA 

The simulations discussed in this paper have been 
made using the CORSIKA program in version 6.000. All 
bugs found in CORSIKA until version 6.014 have been 
corrected also in the extended version. 
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FIG. 5: The density distribution for atmospheric conditions 
of different climatic zones and seasons, plotted as mass overlay 
for a given altitude. In order to enhance the differences, the 
mass overlay for the altitude h is multiplied by exp(/t/8km). 
The differences for altitudes between 20 and 80 km are most 
important for atmospheric particle fluxes, while the differ- 
ences over 80 km are artificially introduced by constraining 
all models to have zero density at 112.83 km (the starting al- 
titude of CORSIKA). For the air shower development this is 
negligible, because the mass overlay at 80 km amounts to less 
than 10" 2 g/cm 2 . 



The simulation of atmospheric particle fluxes with 
CORSIKA starts by selecting the type of primary par- 
ticles and the ranges for the primary energy, zenith and 
azimuth angles, and by fixing the geographical location 
on Earth. The primary energies vary for all simulations 
reported in this paper between the minimum geomag- 
netic cut-off and 10 15 eV. 

The standard CORSIKA version makes use of a planar 
atmospheric model. This is a good approximation as 
long the zenith angle 9 of the particles does not exceed 
70°. The planar atmosphere approximation is used in 
this paper for the calculation of vertical muon fluxes, 
because the experiments are usually limited to muons 
having zenith angles less than 30°. 

For the simulation of the East- West-effect of atmo- 
spheric muons and for all simulations of atmospheric neu- 
trinos the zenith angles must be varied over the complete 
range. These simulations have been made with the so- 
called "curved" version of CORSIKA. Here the curvature 
of the Earth's atmosphere is approximated by sliding and 
tilting plane atmospheres. Each time the horizontal dis- 
placement of a particle exceeds a limit of 6 to 20 km 
(dependent on the altitude), a transition to a new local 
plane atmosphere is performed |o9j . 
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The different primary particles, i.e. protons and he- 
lium nuclei are simulated in separate runs and the ratio 
between them follows the absolute fluxes reported by the 
AMS prototype mission. In order to account for heav- 
ier primary particles the equivalent number of primary 
helium nuclei is used. The absolute fluxes of heavier nu- 
clei are taken from the compilation of Wiebel-Sooth et 
al. |70| . A justification of this simplification is provided 
by the fact that all heavier particles contribute together 
less than 5 % to the neutrino flux and all nuclei have a 
similar ratio of protons to neutrons. 

The air shower calculation starts by getting a random 
location on the Earth, a random energy and a random 
arrival direction. If the particle does not exceed the 
geomagnetic cut-off for the given location or the solar 
modulation, a new set of geographic coordinates, energy 
and arrival angles is used. If the particle fulfills the re- 
quirements, the geomagnetic parameters, the altitude, 
and the atmosphere, are set according to the geograph- 
ical position. Due to the long measuring time of Super- 
Kamiokande, atmospheric models for summer and winter 
are used in equal parts. 

The primary particle is tracked to the first interaction 
point, given by the cross-section of the particle with air. 
The nuclear reaction is handled by the selected hadronic 
interaction model and all secondary particles are tracked 
up to their decay or further interactions. 

The obtained numbers of atmospheric particles have 
to be normalized to the fluxes of primary particles. For 
sake of simplicity the number of primary particles with 
an energy larger than 1000 GeV in the simulation, being 
free of any influence of the geomagnetic cut-off and the 
solar modulation, is set equal to the integral flux above 
1000 GeV as extrapolated in Sec. Illljl In cases with a 
limited statistical accuracy the calibration is made at 
100 GeV. The fluxes at this energy arc already influenced 
by the solar modulation by some 4.5%, what has to be 
taken into account. 

Due to the flat or partially flat geometry applied in 
CORSIKA, the obtained neutrino fluxes have to be scaled 
by the surface difference of the two shells having the 
radius of the Earth and the radius of the Earth plus 
112.83 km. This correction leads to a factor of 1.036. 



TABLE II: The geographical parameters for the different 
detector sites. The quantity h is the elevation over sea-level, 
R c the mean vertical geomagnetic cut-off, B x the horizontal 
component of the magnetic field, B z the vertical downward 
component and B a is the angle between the magnetic and 
geographic north direction. The parameters of the magnetic 
field are valid for the year 2000 and an altitude of 56.4 km. 



site 


h [m] 


Rc [GV] 


B x [ixT] 


Bz [//T] 


B a [deg] 


Bucharest 


85. 


5.6 


21.98 


40.96 


3.64 


Fort Sumner 


1270. 


4.2 


22.89 


44.33 


9.44 


Lynn Lake 


360. 


0.5 


9.81 


57.51 


9.36 


Okayama 


5.3 


11.8 


30.48 


34.31 


-6.64 


Tsukuba 


30. 


11.5 


29.08 


34.82 


-6.95 



comparisons of this communication are focused on the re- 
cent measurements of BESS, CAPRICE, the OKAYAMA 
cosmic ray telescope and WILLI. In case of BESS [7j| 
and CAPRICE [76j the results of atmospheric muons 
have been obtained in ground based runs, performed as 
test of the detectors. The OKAYAMA telescope [71 is a 
classical magnetic spectrometer and WILLI |7l|, [78| rep- 
resents a compact scintillator experiment dedicated to 
the precise measurement of the muon charge ratio. The 
charge ratio is deduced hereby from the different life time 
of positive and negative muons in matter. 

For the simulation of the atmospheric muon fluxes the 
precise geographical parameters, like the geomagnetic 
cut-off and the altitude of the different detector sites are 
taken into account. The used parameters are compiled 
in Tab. [H] Due to the geographic vicinity of Okayama 
and Tsukuba and the same altitude of both sites, the 
results of the OKAYAMA telescope can be compared di- 
rectly with the measurements and calculations done for 
Tsukuba. 

The results for the differential flux of vertical muons 
are compiled in Fig. The calculation with DPMJET 
as well as the calculations with VENUS + UrQMD agree 
generally well with the experimental data. Only the 
GHEISHA results show a strange enhancement of the 
differential muon flux for low energies and quite a differ- 
ent momentum dependence. 



III. CALCULATION OF ATMOSPHERIC MUON 
FLUXES 

A. The differential muon flux 

The calculation of atmospheric muon fluxes controls 
the calculations of atmospheric neutrino fluxes. The 
charge ratio of muons provides additional and partly 
complementary information. 

Atmospheric muons have been measured over sev- 
eral decades. The data are compiled in two recent pa- 
pers 0, [72 and in the new review [7i| , showing rela- 
tively large discrepancies between the experiments. The 



B. The charge ratio of muons 

In contrast to the differential muon fluxes, the charge 
ratio of muons reveals larger discrepancies. The COR- 
SIKA results for the charge ratio of muons are compared 
in Fig. [7| with the experimental data. Again the results 
obtained with the GHEISHA model are far from the ex- 
perimental observations but there are also differences be- 
tween the results of DPMJET and VENUS + UrQMD. 
The results obtained with VENUS + UrQMD are lower 
than the experimental values especially for low and inter- 
mediate energies. It has been shown, that this deviation 
originates mainly from UrQMD while at higher energies 



9 



VENUS leads to a muon charge ratio which is compatible 
to the measurements [79j . 

The DPMJET results agree generally well with the 
data, with exception of the CAPRICE results in Fort 
Sumner. The deviation for Fort Sumner has to be ques- 
tioned because the geomagnetic cut-off in Fort Sumner 
resembles that in Bucharest. Therefore the differences 
in the experimental values and the continuous increase 
of the charge ratio in the CAPRICE measurement for 
Fort Sumner far beyond the geomagnetic cut-off seem to 
indicate experimental problems in this particular mea- 
surement. 

The real influence of the geomagnetic cut-off on the 
muon charge ratio can be seen when comparing the 
CAPRICE and BESS results for Lynn Lake, the WILLI 
results for Bucharest, and the BESS results for Tsukuba. 
At higher energies the ratio stays nearly constant, how- 
ever it decreases, when the geomagnetic cut-off clips the 
high excess of low energy primary protons, as can be ob- 
served in the results for Bucharest and Tsukuba. This 
effect is nicely reproduced by CORSIKA using DPMJET 
as interaction model, while using UrQMD the effect is 



DPMJET 
VENUS + UrQMD 



VENUS + GHEISHA 



^ 40 
jn 

> 30 

'E 

a. 

• 20 

cT 

■a 

3 10 
< 

% 40 



z 

~° 30 



20 
10 
40 
30 
20 
10 



Lynn Lake 




— DPMJET 


J* 




v --- VENUS+UrQMD 
^rv VENUS+GHEISHA 


- m/ 

i i i i i i i_i 


1 [ 


i i i i i— ij i ■ 1 — 


Fort Sumner • 
• / 

• 


•• _ 

, *\ 
















i i i i i i i_l 


i l 


1 1 1 1 1— ij 1 ' 1 - 


: Tsukuba 




• CAPRICE 


/ 


■ 


u T BESS 

■ OKAYAMA 

V ■ 


i 






1 




10 



P, [GeV/c] 



FIG. 6: The differential flux of vertical muons calculated by 
CORSIKA using different models for the description of the 
hadronic interaction in comparison with experimental results 
for various detector sites. In order to enhance the differences 
in the region of interest, the fluxes are multiplied with the 
muon momentum p M . 
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FIG. 7: The charge ratio of muons calculated by CORSIKA 
using different models for the description of the hadronic in- 
teraction compared to experimental data from various detec- 
tors. 



covered by intrinsic problems of the model. 

The systematics of the geomagnetic cut-off shows again 
the problem of the CAPRICE results for Fort Sumner. 
The CAPRICE results have practically the same depen- 
dence on the momentum as the BESS results in Tsukuba, 
where the geomagnetic cut-off is nearly 3 times higher. 

It could be argued that Fort Sumner has an altitude 
of 1230 m above sea level and there could be a strong 
dependence of the charge ratio on the altitude, but the 
CORSIKA simulations include the precise altitude and 
the recent results from BESS show only a weak depen- 
dence of the charge ratio on the altitude. The BESS 
data indicate a 3 % difference between Tsukuba and Mt. 
Norikura which has an altitude of 2770m |80|. 

The disability of GHEISHA, the standard hadronic in- 
teraction model in the detector simulation tool GEANT 
3.21, in reproducing the data of atmospheric muons 
surprises. But in fact serious deficits of GHEISHA 
have already been proved in direct model tests. In 
Refs. [zl M, H2 it has been reported that GHEISHA 
violates the energy, momentum, charge and baryon num- 
ber conservation in the single hadronic interaction. 

At least the energy conservation is also violated on av- 
erage as can be shown by the simulation of extensive air 
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showers with standard CORSIKA. CORSIKA allows to 
summarize all the energy deposited in the atmosphere 
during the shower development. Using GHEISHA as low 
energy hadronic interaction model, an augmentation of 
energy of a complete shower is observed. This increase 
of energy is about 5 % at 10 15 eV and 7 % at 10 14 eV. 
Therefore, the GHEISHA version used in GEANT3 
should not be used in any serious simulation of atmo- 
spheric neutrino fluxes. This holds especially for the 
neutrino flux calculations of Plyaskin, which are based 
on GHEISHA, only. After finishing the simulations, cor- 
rection patches for GHEISHA became available which 
improve essentially the energy conservation |83|. 
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FIG. 8: The East- West-effect of the muon charge ratio as 
measured by the WILLI detector in comparison with calcula- 
tions of CORSIKA using DPMJET. The detector acceptance 
of WILLI is taken into account by processing the raw results 
of CORSIKA by the detector simulation program of WILLI. 
East and WEST mean that the detector looks to the East 
and West, respectively. 



C. The East-West-effect of the muon charge ratio 

Data for inclined muons allow a check of the calcula- 
tions in the curved geometry of the Earth. Using the so- 
called East-West-effect of the muon charge ratio, caused 
by the influence of the geomagnetic field, the way of han- 
dling the field in the calculation can be verified, too. 

Fig. [5] shows preliminary results of the WILLI experi- 
ment for muons observed in East and West direction hav- 
ing a mean zenith angle of 35° 84] in comparison with 
CORSIKA simulations on basis of DPMJET. The COR- 
SIKA results are processed by a full detector simulation 
of the experiment in order to account for the complex 
acceptance of the instrument. 

The agreement of the CORSIKA results with the 
strong East- West-effect observed by the WILLI experi- 
ment, gives confidence that the corresponding effect in 
the atmospheric neutrino flux is also handled well by 
CORSIKA. 



Muon data in various depths of the atmosphere would 
provide a further possibility for the revision of calcula- 
tions on atmospheric particle fluxes. Unfortunately the 
rise and decline time of the actual balloon measurements 
are such fast, that the corresponding muon data have 
large statistical errors. Additionally the atmospheric 
pion flux causes systematic errors in some instruments. 
While the pion flux on sea-level is only 0.5 % of the muon 
flux it reaches 50 % when approaching the top of atmo- 
sphere. 

Nevertheless it has to be pointed out that the atmo- 
spheric muon flux, in contrast to the neutrino flux where 
every produced neutrino reaches ground level, is a highly 
differential quantity, because most muons are already ab- 
sorbed before reaching ground level. Therefore possible 
differences, for example in the nuclear interaction models 
are enhanced from one hadronic interaction to the next. 
Thus the calculation of the ground level muon flux has 
higher theoretical uncertainties than the calculation of 
atmospheric neutrino fluxes. 



IV. CALCULATION OF ATMOSPHERIC 
NEUTRINO FLUXES 

A. The vertical neutrino fluxes in Kamioka 

The calculation of atmospheric neutrino fluxes for 
Kamioka is split in two separate calculations. The down- 
ward going neutrinos are simulated locally for Kamioka, 
while the upward going neutrinos are calculated from pri- 
mary particles distributed over the entire Earth and only 
neutrinos passing in a circle of 1000 km distance from 
Kamioka are used in the further analysis. 
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FIG. 9: The vertical differential intensities of the different 
neutrino flavors in Kamioka, displayed as the ratio between 
the CORSIKA results using DPMJET as hadronic interaction 
model and the calculations of BGS, HKHM and BFLMSR. 
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TABLE III: The vertical differential intensity f(v) and the errors A(V) for downward going neutrinos in Kamioka as calculated 
with CORSIKA using DPMJET and VENUS + UrQMD. The fluxes and the corresponding statistical errors due to the limited 
number of events in the Monte Carlo simulation are given in units of (m 2 s sr GeV)" 1 . 
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This procedure causes a large difference in the number 
of primary particles needed in the simulation for obtain- 
ing the same statistical accuracy for the up- and down- 
ward going fluxes. In the present simulation the number 
of upward going neutrinos is still a factor of 8 smaller. 

Tab. IIIII gives the differential intensities for vertical 
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FIG. 10: The vertical differential intensities of the differ- 
ent neutrino flavors in Kamioka. Shown is the ratio between 
the CORSIKA results using VENUS + UrQMD as hadronic 
interaction model and the calculations of BGS, HKHM and 
BFLMSR. 



neutrinos obtained with CORSIKA, using DPMJET and 
VENUS + UrQMD, respectively. In Figs.E]and[in]the re- 
sults are compared directly with the calculations of BGS, 
HKHM and BFLMSR. 

The inclusive neutrino fluxes obtained with CORSIKA 
are evidently lower than the fluxes given by BGS and 
HKHM. The differential fluxes at 0.1 GeV are about 40 % 
smaller than the BGS fluxes and become comparable 
at energies in the GeV range. The agreement of the 
CORSIKA results using DPMJET and using VENUS + 
UrQMD with the BFLMSR calculation is better. The 
deviation of these absolute flux calculations over prac- 
tically the whole energy range remains less than 20%. 
The energy dependence of the neutrino fluxes between 
BFLMSR and VENUS + UrQMD is quite similar while 
DPMJET shows a systematic difference to BFLMSR. 

Fig. ^3 displays the ratio between the different neu- 
trino flavors in the vertical downward going flux. The 
agreement of all calculations for the ratio of muon neu- 
trinos to electron neutrinos is very good. The deviation of 
the HKHM results and the discontinuity at E v = 1 GeV 
are caused by different approaches in the model. Below 
1 GeV the values of HKHM are averaged over the zenith 
angle, only above 1 GeV they stand for vertical, down- 
ward going neutrinos. For energies below 3 GeV the dif- 
ferences between the other models are on the level of 2 % 
or better. 

Some differences between the calculations are observed 
in the ratio of muon neutrinos to muon antineutrinos. 
The CORSIKA calculations with DPMJET and VENUS 
+ UrQMD, and the BFLMSR calculations agree per- 
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fectly. The calculations ol BGS predict a lower ratio 
at 3 GeV while the calculations of HKHM are different 
around 1 GeV and show a smaller rise of the ratio at high 
energies. 

The ratio of electron neutrinos to electron antineutri- 
nos reveals larger differences. The results of HKHM be- 
have quite different from the results of all other models. 
Interestingly, DPMJET results agree with BFLMSR re- 
sults, while VENUS + UrQMD results agree with BGS 
results. Due to the close correlation between the ratio 
of electron neutrinos to electron antineutrinos and the 
charge ratio of muons, these findings allow to rule out the 
results of VENUS + UrQMD in this particular quantity, 
meaning that the results of BGS are suspicious in this 
aspect, too. 
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FIG. 11: The ratio between different neutrino flavors in the 
vertical flux in Kamioka as calculated by CORSIKA with 
DPMJET and with VENUS + UrQMD. 

An interesting quest for the CORSIKA calculations 
with their inclusion of the precise geometry of the Earth, 
are natural differences between the up- and downward 
going neutrino fluxes in Kamioka. Such differences could 
contribute to the measured asymmetry, which is com- 



monly attributed to the oscillation of neutrinos. Any 
natural difference based on the geographical environment 
has a direct impact on the analysis of the neutrino oscil- 
lations and changes finally the obtained oscillation pa- 
rameters. 
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FIG. 12: The asymmetry of the up- (N u ) and downward (Nd) 
going neutrino fluxes for Kamioka as calculated by CORSIKA 
using DPMJET and VENUS + UrQMD in comparison with 
the calculation of BFLMSR. The asymmetry is expressed by 
the ratio between the difference and the sum of the up- and 
downward going neutrino fluxes. 

A major difference between Kamioka and its antipode 
in the South Atlantic comes from the geomagnetic cut- 
off. While the vertical cut-off in Kamioka is 12.3 GV, 
the South Atlantic region is influenced by the so-called 
South Atlantic magnetic field anomaly leading to a verti- 
cal cut-off at the antipode of only 8.6 GV. This causes an 
asymmetry between the intensities of up- and downward 
going neutrinos for Kamioka, as can be seen in Fig. ^| 

The asymmetry of 20 % having been observed in the 
calculations of BFLMSR represents the raw effect based 
on the differences in the geomagnetic cut-off, because the 
calculation does not include any local magnetic field. In 
the CORSIKA simulations the local field and an addi- 
tional contribution to the up-down asymmetry, caused 
by the different elevation of the surface over sea-level in 
Kamioka and in the antipode region in the South At- 
lantic, are taken into account. 

The location of the Super-Kamiokande detector in the 
mountains causes an altitude difference of several hun- 
dred meters compared to the average altitude of the an- 
tipode region. Thus in the South Atlantic the shower 
development is longer and more neutrinos are produced 
in the shower. Further details on the influence of the 
local magnetic field and the geomagnetic cut-off are in- 
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FIG. 13: The zenith angle dependence of the neutrino intensities in Kamioka as calculated by CORSIKA with DPMJET and 
with VENUS + UrQMD in comparison with the calculations of BFLMSR and BGS. The first column shows the plots for the 
energy interval 0.1 to 0.31 GeV, the second for 0.31 to 1 GeV, the third for 1 to 3.1 GeV, and the last for 3.1 to 10 GeV. 



vestigated in Sec. lIV Cl The effect of the contrary seasons 
in Japan and the South Atlantic, which is taken into ac- 
count by using the appropriate atmospheric models does 
not lead to any observable effect, the effect is smaller 
than the actual statistical errors. 



B. The directional dependence of the neutrino 
fluxes in Kamioka 

The dependence of the neutrino fluxes on the zenith 
angle is shown in Fig. The three dimensional calcula- 
tions of BFLMSR and CORSIKA show an enhancement 
of the neutrino fluxes near the horizon. This enhance- 
ment is based on a geometrical effect, i.e. the spherical 
shell geometry of the neutrino production volume |85j . 
This effect has been neglected in all one dimensional sim- 
ulations like HKHM and BGS. The strength of the effect 
shows clearly the necessity of the time consuming three 
dimensional simulations in a spherical geometry. The 
agreement of the calculation with VENUS + UrQMD 



and with BFLMSR is again better, while the DPMJET 
results show systematically higher fluxes for energies be- 
tween 1 and 3 GeV. 

The dependence of the resulting ratio between muon 
neutrinos and electron neutrinos on the zenith angle is 
shown in Fig. Only the results for energies below 
1 GeV arc plotted, for higher energies no difference be- 
tween all the four calculations is observed. As in the case 
of the ratios between vertical neutrino fluxes the largest 
differences are observed in the ratio of electron neutrinos 
to electron antineutrinos. The CORSIKA results show a 
strong increase of the ratio near the horizon. The origin 
of this effect will be investigated in Sec. II V CI 

Also a 8 % difference of the ratio of muon neutrinos to 
electron neutrinos at low energies can be observed near 
the horizon. The results of the BGS calculation lead to 
very low values for this quantity, and may be an artifact 
of the calculation in a one dimensional geometry. 

The dependence of the neutrino fluxes on the azimuth 
angle is shown in Fig. 1151 The agreement between the 
calculations with DPMJET and with VENUS + UrQMD 
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FIG. 14: The results of CORSIKA using DPMJET and 
VENUS + UrQMD for the zenith angle dependence of the 
different ratios between the neutrino flavors. The values are 
compared with the calculations of BFLMSR and BGS. The 
plots in the first row are for the energy interval 0.1 to 0.31 GeV 
and in the second row for 0.31 to 1 GeV. If the error bar is 
not drawn, the error is smaller than the symbol size. 

for westward going neutrinos is very good, but for east- 
ward going neutrinos some noticeable differences are ob- 
served at higher energies. This is a secondary effect of 
the difference in the momentum spectra of the reaction 
products between both models, but it displays also an 
instructive example how the interaction model influences 
results which are commonly assumed to have a geomet- 
rical nature. 

The detailed comparison with the results of the HKKM 
calculation shows a very good agreement in the shape 
of the azimuthal distribution. At lowest energies the 
HKKM calculation leads to much higher fluxes. The au- 
thors state this overestimation to be caused by the use 
of the old COSMOS interaction models. A new calcula- 
tion using DPMJET as hadronic interaction model will 
overcome this problem. 

The good agreement between CORSIKA results and 
the calculation of HKKM in the azimuthal distribution 
is by far not trivial, as shows the comparison of the COR- 
SIKA results with calculations of Lipari et al. [85ll86| in 
Fig. ^3 Here the shapes of the distributions for electron 
neutrinos and muon antincutrinos are compatible, but 
strong disagreement exists for electron antineutrinos and 
muon neutrinos. 

The results can be expressed by the East- West- 
asymmetry Aew = (Ne — Nw)/{Ne + Nw), where Ne 
and Nyy stand for the particle fluxes of neutrinos going 
to the East and West, respectively. Fig. El shows the 
energy dependence of the East- West-asymmetry. Again 



the CORSIKA results with DPMJET have a slightly 
higher asymmetry than the calculations with VENUS + 
UrQMD. The distributions of all neutrino flavors show 
similar shapes. The strongest asymmetry is observed for 
electron neutrinos and the weakest for electron antincu- 
trinos. All neutrino flavors exhibit a maximal asymmetry 
for an energy around 800 MeV. 



C. The influences of the geomagnetic cut-off and 
the local magnetic field 

In order to investigate the individual influences of 
the geomagnetic cut-off and of the local magnetic field, 
the calculations of the atmospheric neutrino fluxes for 
Kamioka with DPMJET have been repeated twice under 
the same conditions, except setting once the local mag- 
netic field and once the geomagnetic cut-off to zero. This 
procedure allows to disentangle the individual influences 
of the two effects. 

Due to the fact that charged particles do not win or 
loose energy in a magnetic field, the influence of the local 
magnetic field on the total neutrino fluxes is negligible. 
The main effects are expected in the ratios of neutrinos 
and in the azimuthal distribution of the fluxes. Especially 
the ratio of electron neutrinos to electron antineutrinos 
shows a strong effect because the electron neutrinos are 
predominantly produced by positive muons and the elec- 
tron antineutrinos by negative muons. 

Muon neutrinos and muon antineutrinos are produced 
also in the decay of charged pions. fn contrast to the 
muon decay, muon neutrinos result here from the decay 
of positive and muon antineutrinos from the decay of 
negative particles. Due to the shorter life time and the 
higher momentum the total bending of pions is less and 
the bending of the muons is preponderating, but the total 
effect of the local magnetic field on the muon neutrinos 
remains weaker. 

The effect of the inclusion of the local magnetic field in 
the calculation is shown in Fig. ED The increase of the 
ratio between electron neutrinos and electron antineutri- 
nos near the horizon as observed in Fig. ^] has to be at- 
tributed completely to the bending of the charged shower 
particles in the atmosphere. 

The CORSIKA results for the azimuthal dependence of 
the atmospheric neutrino fluxes under the different con- 
ditions are displayed in Fig. 1191 The differences are pro- 
nounced for smaller energies. At higher energies all the 
different conditions lead to identical fluxes. The influence 
on the shape of the azimuthal distribution is weak, but 
only for detector sites with a high geomagnetic cut-off. 

Without consideration of the geomagnetic cut-off, 
much higher neutrino fluxes are obtained due to the 
higher fluxes of primary particles. The asymmetry in 
the azimuthal distribution results here only from the de- 
flection of charged shower particles in the local magnetic 
field. The characteristics of this asymmetry is very sim- 
ilar to the East- West-effect caused by the geomagnetic 
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FIG. 15: The azimuth angle dependence of the neutrino fluxes in Kamioka as calculated with CORSIKA for different energies 
and neutrino flavors compared with the calculations of HKKM. The data in the first row are integrated in an energy interval 
from 0.1 to 0.3 GeV, in the second row from 0.3 to 1 GeV, and in the third row from 1 to 3.1 GeV. The neutrinos selected 
are from both hemispheres and have a | cos#| > 0.5. In this diagram cj> = indicates a particle going in the magnetic North 
direction. The errors are smaller than the symbol sizes. 



cut-off, a consequence of the excess of positive particles 
in the atmosphere, on which the magnetic field acts in 
a similar way as on the primary proton flux. This argu- 
ment is supported by the different behavior of electron 
antineutrinos, which are produced only in the decay of 
negative muons. 

In order to illustrate the transition between a zero and 
a high geomagnetic cut-off, the results of a calculation as- 
suming an isotropic cut-off of 6 GV have been added also 
in Fig. El These results show that a neglect of the local 
magnetic field, as it is done in many calculations of at- 
mospheric neutrino fluxes, may lead to wrong azimuthal 
distributions at least for detector sites with a comparable 
low geomagnetic cut-off. 



V. CONCLUSION 

This work aims at a new procedure for the calculation 
of atmospheric neutrino fluxes with considering various 
influences which have not been taken into account so far 
or, if ever, only in a less rigorous way. The capabilities 
of the procedure are demonstrated by a particular cal- 
culation of the detailed neutrino fluxes in Kamioka. The 
detailed procedure applies the air shower simulation code 
CORSIKA in the version 6.000, which has been extended 
and modified for a reliable simulation of the cascading 
interactions induced in the atmosphere by primary par- 
ticles of the low energy part of the cosmic ray spectrum. 

A description of the solar modulation and tables for 
the geomagnetic cut-off, calculated in a detailed Monte 
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FIG. 17: The energy dependence of the East-West- 
asymmetry in the atmospheric neutrino flux as calculated by 
CORSIKA with DPMJET and with VENUS + UrQMD. 



Carlo simulation, have been introduced. In addition, for 
the first time for atmospheric neutrino flux calculations, 
the geography of the Earth is taken into account by a 
digital elevation model, tables for the local magnetic field 
in the atmosphere, and various atmospheric models for 
different climatic zones and seasons. These extensions 
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FIG. 18: Zenith angle dependence of the ratio between elec- 
tron neutrino and electron antineutrino fluxes in Kamioka. 
The results of CORSIKA with DPMJET considering the lo- 
cal magnetic field are compared with results of CORSIKA 
and BFLMSR neglecting the local magnetic field. 



are not yet part of the standard CORSIKA package. 

CORSIKA features a precise particle tracking, includ- 
ing the deflection of the charged shower particles in the 
local magnetic field, the energy loss by ionization and 
multiple scattering. The used primary flux is based on 
the recent measurements of the prototype of the AMS- 
experiment. These data allow also a test of the calcula- 
tions for the geomagnetic cut-off. 

An important aspect of the calculations is the question 
of the adequate hadronic interaction model used as gener- 
ator of the flux calculations. This question is approached 
by using the possibilities of the CORSIKA code to op- 
erate optionally with various different hadronic interac- 
tion models. The models are scrutinized by an extensive 
comparison with measured fluxes and charge ratios of at- 
mospheric muons in different locations. 

It turns out that the GHEISHA model leads to signif- 
icant discrepancies with data from various experiments 
and predictions based on GHEISHA have to be consid- 
ered as highly doubtful. The use of DPMJET II.5 as well 
as of the combination VENUS + UrQMD results in dif- 
ferential muon fluxes which are in good agreement with 
the measurements. The DPMJET model reproduces the 
charge ratio of muons of vertical incidence, while the val- 
ues obtained with VENUS + UrQMD appear systemat- 
ically too small. The calculations with DPMJET agree 
also well with the preliminary results of the WILLI ex- 
periment for the East- West-effect of the charge ratio of 
muons with inclined incidence. 

Subsequently CORSIKA is used with the described re- 
finements to calculate the fluxes of atmospheric neutrinos 
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FIG. 19: The azimuth angle dependence of the atmospheric neutrino fluxes in Kamioka as calculated by CORSIKA under 
different conditions. The neutrinos used in the analysis result from both hemispheres and have a |cos#| > 0.5. Beside the 
results described already above, the azimuthal distributions are shown for a zero local magnetic field, a zero geomagnetic cut-off 
and an isotropic cut-off of 6 GV. In this diagram <j> = indicates a particle going towards the magnetic North direction. The 
distributions in the first row are integrated over an energy interval between 0.1 and 0.3 GeV, in the second between 0.3 and 
1 GeV, and in the third between 1 and 3.1 GeV. 



for Kamioka. The resulting absolute neutrino intensities 
are lower than those found in the classical calculations of 
BGS and HKHM, but they are in good agreement with 
the recent three dimensional calculations of BFLMSR. 
Using VENUS + UrQMD the deviations from BFLMSR 
predictions are smaller than 20 % over the whole energy 
range and the overall energy dependence is very similar. 

DPMJET leads to absolute fluxes, being also very simi- 
lar to the simulations of BFLMSR, but the energy depen- 
dence turns out to be slightly different. Nevertheless the 
better agreement of the DPMJET predictions with the 
measured fluxes and charge ratios of atmospheric muons 
provides stringent arguments in favor of this particular 
model. 

The ratio of muon neutrinos to electron neutrinos and 



the ratio of muon neutrinos to muon antineutrinos in 
the vertical downward flux are identical within the sta- 
tistical uncertainties for the CORSIKA calculations in- 
voking DPMJET, VENUS + UrQMD, and the calcula- 
tions of BFLMSR. But for lowest energy neutrinos with 
horizontal incidence, the ratios between muon neutrinos 
and electron neutrinos obtained with DPMJET and with 
VENUS + UrQMD are higher. 

Significant differences are observed for the ratio of 
electron neutrinos to electron antineutrinos. The DPM- 
JET results for vertical neutrinos for this quantity agree 
with the results of BFLMSR, and the results of VENUS 
+ UrQMD agree with the results of BGS. Again the 
very good agreement in the correlated quantity of the 
muon charge ratio gives a strong argument for DPMJET. 
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For horizontal neutrinos the CORSIKA results predict a 
strong increase of the ratio at low energies. 

The actual results have relevance for the analysis of 
the atmospheric neutrino anomaly. Any change in the 
ratio of muon neutrinos to electron neutrinos leads di- 
rectly to a change of the oscillation parameters. Also 
the discrepancies found in the ratio of electron neutri- 
nos to electron antineutrinos are of particular interest for 
Super-Kamiokande, because the detection cross-sections 
for neutrinos are about three times larger than for an- 
tineutrinos and it is not possible to distinguish between 
them in the experiment. 

To quantify the influence of these effects on the neu- 
trino oscillation parameters would request a full detector 
simulation of the Super-Kamiokande experiment based 
on the presented fluxes, a task which is beyond the scope 
of this communication. It can be stated that the differ- 
ence of the neutrino fluxes presented here to those used 
in the oscillation analysis is not large enough to affect 
the claim of existence of neutrino oscillations from atmo- 
spheric neutrinos. 

The use of two different hadronic interaction models, 
both of good repute in interpretation of accelerator ex- 
periments, shows clearly the potential influence of the 
hadronic interaction model on the interpretation of the 
atmospheric neutrino anomaly. Due to the high quality 
of the recent measurements of the primary particle fluxes, 
the main source of remaining uncertainties in the atmo- 
spheric flux calculations has to be attributed now to the 
actual uncertainties in the hadronic interaction models. 

For studying the influence of the geomagnetic field and 
the origin of the East- West-effect in the atmospheric neu- 
trino flux, CORSIKA calculations with DPMJET, set- 
ting the local magnetic field to zero or skipping the ge- 
omagnetic cut-off have been performed. The main influ- 
ence of the local magnetic field is found for the ratio of 
electron neutrinos to electron antineutrinos. CORSIKA 
predicts for the first time a strong increase of the ratio 
near the horizon. 

The local magnetic field proves to be of minor in- 
fluence on the azimuthal distribution of neutrinos in 
Kamioka, and the East- West-effect arises mainly from 
the azimuthal dependence of the primary particle flux 
caused by the geomagnetic cut-off rigidity. The simula- 



tions without a geomagnetic cut-off show that this obser- 
vation is valid only for Kamioka with its relative high ge- 
omagnetic cut-off value. For a neutrino detector site like 
Sudbury in Canada, where the vertical geomagnetic cut- 
off is only 1.1 GV, a measurable East- West-effect would 
originate exclusively from the bending of the charged 
shower particles in the local magnetic field. 

To which extent the Earth's geography significantly 
affects the results of the calculations has not been inves- 
tigated in detail by separate calculations. The higher 
asymmetry of the up- and downward going particle 
fluxes, found in the actual calculations in comparison to 
results of BFLMSR, indicates an influence of the digital 
elevation model in the order of a few percent. Compared 
to the changes of the atmospheric depth by the different 
altitudes, the variation induced by the different atmo- 
spheric models is small. The influence on the particle 
fluxes in Kamioka should be negligible. Nevertheless for 
detector sites with extreme atmospheric conditions, like 
the South Pole, the profile of the atmosphere may lead 
to noticeable seasonal effects. 



Acknowledgments 

This work has been supported by the Deutsche 
Forschungsgcmeinschaft with the grants WE 2426/ Tl 
and WE 2426/1-2. The help of the International Bu- 
reau Bonn supporting the personal exchange and of the 
Volkswagen-Stiftung for sponsoring valuable devices is 
gratefully acknowledged. The authors are deeply in- 
debted to G. Schatz and H. Bliimer which enabled 
and supported the major part of this study in the 
Forschungszcntrum Karlsruhe. The suggestions and 
valuable advice of H. Stocker and M. Bleicher when in- 
corporating the UrQMD model in the CORSIKA code 
and of J. Ranft when applying DPMJET at low energies, 
are appreciated. We thank R. Engel for carefully reading 
the manuscript and T.K. Gaisser for providing us with 
tables of the BGS results. One of the authors (JW.) 
is grateful to the European Commission Centre of Ex- 
cellence (IDRANAP) in Bucharest for the grant which 
allowed the finalizing of this paper. 



[1] Y. Fukuda et al., Phys. Lett. B 433, 9 (1998). 
[2] Y. Fukuda et al., Phys. Rev. Lett. 81, 1562 (1998). 
[3] R. Becker-Szendy et al., Phys. Rev. D 46, 3720 (1992). 
[4] K. S. Hirata et al., Phys. Lett. B 280, 146 (1992). 
[5] W. W. M. Allison et al., Phys. Lett. B 391, 491 (1997). 
[6] T. Kafka, Nucl. Phys. (Proc. Suppl.) B 70, 340 (1999). 
[7] T. Futagami et al., Phys. Rev. Lett. 82, 5194 (1999). 
[8] IAGA Division V, Working Group 8, AOS, Trans. AGU 

73, 182 (1992). 
[9] S. MacMillan et al., Jour, of Geomag. and Geoelec. 49, 

229 (1997). 



[10] J. M. Quinn et al., Jour, of Geomag. and Geoelec. 49, 
245 (1997). 

[11] S. E. Forsythe, Smithsonian Physical Tables, Smithso- 
nian Institution Press (1969). 

[12] T. K. Gaisser, T. Stanev, and G. Barr, Phys. Rev. D 38, 
85 (1988). 

[13] G. Barr, T. K. Gaisser, and T. Stanev, Phys. Rev. D 39, 
3532 (1989). 

[14] V. Agrawal et al., Phys. Rev. D 53, 1314 (1996). 
[15] E. V. Bugaev and V. A. Naumov, Phys. Lett. B 232, 391 
(1989). 



19 



Phys. Lett. B 248, 193 (1990). 
Phys. Rev. D 52, 4985 (1995). 
Bludman, Phys. Rev. D 37, 122 (1988). 
Koh, Nuovo Cim. B 105, 883 (1990). 

(1999), hep-ph/9907450; As- 

yth 



M. Honda et al., 
M. Honda et al., 
H. Lee and S. A. 
H. Lee and Y. S. 
Y. Tserkovnyak et al. 
tropart. Phys., in press. 

Y. Tserkovnyak et al., in Proc. 27 th Int. Cosmic Ray 
Conf. (Hamburg, 2001), p. 1196. 

G. Battistoni et al., Astropart. Phys. 12, 315 (2000). 
G. Battistoni, Nucl. Phys. B (Proc. Suppl.) 100, 101 
(2001). 

G. Battistoni et al. (2002), hep-ph/0207035. 

M. Honda et al., Phys. Rev. D 64, 053011 (2001). 

M. Honda et al., in Proc. 27 th Int. Cosmic Ray Conf. 

(Hamburg, 2001), p. 1162. 

V. Plyaskin, Phys. Lett. B 516, 213 (2001). 

T. K. Gaisser and M. Honda, Ann. Rev. of Nucl. Part. 

Sci. 52, 153 (2002). 

T. K. Gaisser, R. J. Protheroe, and T. Stanev, in Proc. 
18 th Int. Cosmic Ray Conf. (Bangalore, 1983), vol. 5, p. 
174. 

K. Kasahara, in Proc. 24 th Int. Cosmic Ray Conf. 
(Rome, 1995), vol. 1, p. 399. 

B. Nilsson-Almqvist and E. Stenlund, Comp. Phys. Com- 
mun. 43, 387 (1987). 

T. Sjostrand and M. Bengtsson, Comp. Phys. Commun. 
43, 367 (1987). 

K. Hanfigen and J. Ranft, Comp. Phys. Commun. 39, 37 
(1986). 

GEANT Detector Description and Simulation Tool, 
CERN (1993), Program Library Long Writeup W5013. 

H. W. Bertini, Phys. Rev. 131, 1801 (1963). 
H. W. Bertini, Phys. Rev. C 1, 423 (1970). 
H. W. Bertini, Phys. Rev. C 6, 631 (1972). 

P. A. Aarnio et al., in Proc. of the MC93 Int. Conf. 

on Monte Carlo Simulation in High Energy and Nuclear 

Physics, World Scientific (Tallahassee, 1993), p. 88. 

A. Fasso et al., in Proc. of the Workshop on Simulating 

Accelerator Radiation Environment (Santa Fe, 1993), vol. 

LA-12835-C of Los Alamos Reports, p. 134. 

H. Fesefeldt, Report PITHA-85/02, RWTH Aachen 

(1985). 

A. Fasso et al., in Proc. of the Monte Carlo 2000 Conf., 

Springer Verlag (Lisbon, 2000), p. 159. 

A. Fasso et al., in Proc. of the Monte Carlo 2000 Conf., 

Springer Verlag (Lisbon, 2000), p. 955. 

S. Roesler, R. Engel, and J. Ranft (2000), hep- 

ph/0012252. 

D. Heck et al., Report FZKA 6019 (1998). 

J. Wentz et al., in Proc. 27 th Int. Cosmic Ray Conf. 

(Hamburg, 2001), p. 1167. 

J. Ranft, Phys. Rev. D 51, 64 (1995). 

J. Ranft (1999), hep-ph/9911213. 

J. Ranft (1999), hep-ph/9911232. 

K. Werner, Phys. Rep. 232, 87 (1993). 

N. N. Kalmykov, S. S. Ostapchenko, and A. I. Pavlov, 

Nucl. Phys. B (Proc. Suppl.) 52B, 17 (1997). 

J. Engel et al., Phys. Rev. D 46, 5013 (1992). 

R. S. Fletcher et al., Phys. Rev. D 50, 5710 (1994). 



[53' 
[54 

[55' 
[56' 

[57; 

[58 
[59' 
[60 
[61 

[62; 

[63 

[64 

[65 
[66' 
[67 

[68' 



[69 

[70; 

[71 

[72; 
[73; 

[74 

[75; 

[76 
[77 

[78; 

[79 

[80 
[81 

[82; 

[83; 

[84 

[85' 
[86; 



ig. Part. Nucl. Phys. 
Phys. G: Nucl. Part. 



41, 255 (1998). 
Phys. 25, 1859 



Phys. Rev. D 60, 052002 (1999). 
Astrophys. Jour. 545, 1135 (2000). 
Astrophys. Jour. 518, 457 (1999). 
Astrophys. Jour. 533, 281 (2000). 
Phys. Lett. B 490, 27 (2000). 
Phys. Lett. B 494, 193 (2000). 



S. A. Bass et al., Pr< 
M. Bleicher et al., J. 
(1999). 

R. Belloti et al., 
T. Sanuki et al. 
M. Boezio et al.. 
W. Menn et al., 
J. Alcaraz et al. 
J. Alcaraz et al.. 
T. K. Gaisser et al., in Proc. 27 th Int. Cosmic Ray Conf. 
(Hamburg, 2001), p. 1643. 

D. E. Page and R. G. Marsden, eds., The Heliosphere at 
Solar Minimum (Pergamon Press, 1997). 
L. J. Gleeson and W. I. Axford, Astrophys. Jour. 154, 
1011 (1968). 

J. Wentz, A. Bercuci, and B. Vulpescu, in Proc. 27 th Int. 

Cosmic Ray Conf. (Hamburg, 2001), p. 4213. 

J. Alcaraz et al., Phys. Lett. B 472, 215 (2000). 

P. Lipari, Astropart. Phys. 16, 295 (2002). 

National Geophysical Data Center, Data Announcement 

88-MGG-02 (1988). 

F. X. Kneizys et al., The MODTRAN 2/3 Report and 
LOWTRAN 7 Model, Phillips Laboratory, Hanscom 
AFB, Massachusetts (1996). 

D. Heck et al., in Proc. 26 th Int. Cosmic Ray Conf. (Salt 
Lake City, 1999), vol. 1, p. 498. 

B. Wiebel-Sooth, P. L. Biermann, and H. Meyer, Astron. 
and Astrophys. 330, 389 (1998). 

B. Vulpescu et al., Nucl. Instr. and Meth. A 414, 205 
(1998). 

T. Hebbeker and C. Timmermans, Astropart. Phys 18, 
107 (2002). 

P. K. F. Grieder, Cosmic Rays on Earth, Researcher's 
Reference Manual and Data Book (Elsevier, 2001). 



in Proc. 27 Int. Cosmic Ray Conf. 
p. 927. 

(2002), astro-ph/0205344; Astropart. 



83, 4241 (1999). 
Part. Phys. 24, 



1805 



Part. Phys. 27, 977 
Phys. 27, 1699 



M. Motoki et al. 
(Hamburg, 2001) 
M. Motoki et al. 
Phys., in press. 

J. Kremer et al., Phys. Rev. Lett. 
S. Tsuji et al., J. Phys. G: Nucl 
(1998). 

B. Vulpescu et al., J. Phys. G: Nucl 
(2001). 

J. Wentz et al., J. Phys. G: Nucl. Part 
(2001). 

T. Sanuki et al., Phys. Lett. B 541, 234 (2002). 

A. Ferrari and P. R. Sala (1996), Atlas Internal Note, 

Phys-Np-086, unpublished. 

J. Wentz et al., in Proc. 26 th Int. Cosmic Ray Conf. (Salt 
Lake City, 1999), vol. 2, p. 92. 

R. E. Cassell and G. Bower (2002), private communica- 
tion. 

I. M. Brancus et al., in Proc. of the 16 th Particles and 
Nuclear Int. Conf. (Osaka, 2002), Nucl. Phys. A, in press. 
P. Lipari, Astropart. Phys. 14, 171 (2000). 
P. Lipari, T. Stanev, and T. K. Gaisser, Phys. Rev. D 
58, 073003 (1998). 



